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ABSTRACT 

This paper examines how to compare stellar limb-darkening coefficients evaluated 
from model atmospheres with those derived from photometry. Different characteriza- 
tions of a given model atmosphere can give quite different numerical results (even for 
a given limb-darkening 'law'), while light-curve analyses yield limb-darkening coeffi- 
cients that are dependent on system geometry, and that are not directly comparable 
to any model-atmosphere representation. These issues are examined in the context of 
exoplanetary transits, which offer significant advantages over traditional binary-star 
eclipsing systems in the study of stellar limb darkening. 'Like for like' comparisons be- 
tween light-curve analyses and new model-atmosphere results, mediated by synthetic 
photometry, are conducted for a small sample of stars. Agreement between the re- 
sulting synthctic-photometry/atmosphere-model (SPAM) limb-darkening coefficients 
and empirical values ranges from very good to quite poor, even though the targets 
investigated show only a small dispersion in fundamental stellar parameters. 
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1 INTRODUCTION 

Stellar limb darkening is the wavelength-dependent decrease 
in specific intensity, I\((J,), with decreasing fi, where p = 
cos 9 and 6 is the angle between the surface normal and the 
line of sight[J]in the context of model atmospheres, it is, in 
principle, significantly more sensitive to input physics than 
are integral quantities, such as the emergent flux. 

Until rather recently, the only important opportunity to 
compare models and observations of limb darkening for the 
distant stars has been through eclipsing-binary systems, but 
there the comparison has been hindered both by the rather 
weak dependence on limb darkening of the light-curves, and 
by degeneracies with other model parameters. As a con- 
sequence, normal practice among light-curve analysts has 
been to assume some description of limb darkening, based 
on stellar-atmosphere results; any errors in the description 
are liable to be concealed by small adjustments to fitted free 
parameters. 

New observational techniques have begun to allow the 
direct investigation of limb darkening (and hence more sen- 
sitive tests of model-atmosphere calculations) under other 
circumstances. Optical interferometry has opened the way 
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1 Under some circumstances, limb brightening can occur. 

2 The specific intensity, referred to as the radiance in other con- 
texts, is the rate of energy flow per unit area, per unit time, per 
unit wavelength interval, per unit solid angle. Expressing I\ as 
a function of a single angle fi makes the implicit assumption of 
azimuthal symmetry of the radiation field. 



to direct imaging of stellar surfaces beyond the solar system 
for a handful of stars with the largest angular diameters 
(e.g., |Aufdenberg, Ludwig fc Kerve lla 2005), and microlens- 



ing light-curves are also capable of probing the intensity 



distribution of the lensed source (e.g., Witt 


1995||Zub et al. 


2011 1, albeit usually only crudely (Dominik 


2004}. However, 



the focus of the present paper is on the role of limb darken- 
ing in exoplanetary transits, which are likely to yield many 
more results in the coming years than any other technique. 

In many respects, star+exoplanet systems are close to 
being idealised eclipsing binaries: it is often reasonable to 
assume that the photometric properties of the parent star 
are unaffected by the transiting planet (i.e., no tidal distor- 
tion, 'reflection' effect, or gravity darkening), and that the 
secondary (planet) is completely dark, and spherical. These 
assumptions reduce the number of geometric unknowns to 
be determined from the light-curve to only three (in ad- 
dition to the orbital ephemeris, which may be established 
separately); e.g., the ratio of the radii, the size of the star 
in units of the centres-of-mass separation, and the impact 
parameter. This relative simplicity allows a more critical 
examination of limb darkening than is possible in star+star 
systems. With an anticipated torrent of data of extremely 
high quality from satellites such as Kepler, it is therefore 
timely to revisit the comparison of limb-darkening coeffi- 
cients (LDCs) from model-atmosphere and light-curve anal- 
yses, as has already been recognized by several authors (e.g. 
Southworth|2008l |Pal|2008l |Claret|2009b . 

I f I I f I _l r 

This comparison is examined here as follows: Section [2] 
reviews limb-darkening 'laws' and fitting techniques (in- 
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eluding a new flux-conserving least-squares methodology), 
stressing the spread in numerical coefficients that can arise 
even when characterizing a given model-atmosphere inten- 
sity distribution with a given law. Section [3] examines the 
LDCs extracted from light-curve analyses, emphasizing not 
only the range in numerical coefficients that can arise from 
characterizing a given surface-brightness distribution under 
different geometries, but also that the photometrically de- 
termined LDCs are not, in any case, directly comparable to 
those derived from model-atmosphere calculations. 

With the background that (i) the numerical values of co- 
efficients determined from model atmospheres depend on the 
fitting method, and (ii) coefficients determined from light- 
curves are not directly comparable to model-atmosphere re- 
sulst, and vary with impact parameter, Section [4] outlines 
how model-atmosphere limb-darkening results can be com- 
pared with inferences from high-quality exotransit photom- 
etry, and presents illustrative results. 



2 CHARACTERIZING LIMB DARKENING 

In photometric analyses, it is still impractical to invert 
observed light-curves in order to recover detailed stellar 
surface-brightness distributions. Rather, in this context limb 
darkening is habitually represented by some ad hoc 'law' 
with one or, at most, two free parameters, which may be op- 
timized as part of the fitting process]^] In order to facilitate 
comparison of light-curve results with model-atmosphere 
calculations, model intensities are often represented with the 
same parametrizations. 



2.1 Functional forms 

Historically, the first limb-darkening law to be developed 
was the analytical solution for an atmosphere in which the 
source function is linear in optical depth: 



J A (/i) = J A (l) [1-u(1-/0] 



(1) 



( Schwarzschild 1906), where the wavelength dependence of 
u is implicit (although u = 0.6 at all wavelengths for a grey 
atmosphere; Milne|[l~921 1. This linear law is the universally 
adopted one-parameter representation of limb darkening. 

More-realistic atmosphere models do not have analyti- 
cal functional representations of actual limb darkening. Fol- 
lowing the work of IKopall (119491) , a quadratic law of the form 



7(/x) = 7(l) [1-mi(1-/x)- W2 (1- M ) 2 ] 



(2) 



has been widely adopted as a characterization of model- 
atmosphere calculations. It is of particular importance 
in modelling exotransit photometry using Monte-Carlo 
Markov-Chain (MCMC) techniques, since it allows for ana- 
lytical calculation of light-curves with good computational 
efficiency ( jMandel fe Agol | 2002| . 

While eqtns. |1| and |2j are convenient in the analysis 



of light-curves]^] a significantly more accurate representa- 
tion of model-atmosphere results is achieved with the four- 



coefficient fit introduced by Claret (20001 



j(aO = /(i) 



n/2 



(3) 



This form reproduces intensities from model atmospheres 
to ~1 part in 1000 over a wide range of parameter space 
(e.g., Howarth]|2011 1, although it isn't practical to estimate 



numerical values of the coefficients from photometry. 



2.2 Fitting model-atmosphere intensities. 

Although linear and quadratic limb-darkening laws may not 
give particularly accurate functional descriptions of model- 
atmosphere intensities, it is nonetheless necessary to repre- 
sent them in this way in order to compare with observation- 
ally derived LDCs. However, even for a given limb-darkening 
law, the characterization of model-atmosphere results using 
different fitting techniques can result in quite different values 
for the coefficients 



2.2.1 LSI: least squares with 7(1) constrained 
Rewriting eqtn. [T] as 

7 a ( M )//a(1) = [1-u(1-m)], 

gives a one-parameter formulation straightforwardly solved 
by least squares for u, using as input the model-atmosphere 
values of 7(/Lt). The intercept of the linear fit is implicitly con- 
strained such that 7(1), the value of 7(1) evaluated from the 
fitted law, is fixed at model-atmosphere value. The quadratic 
equivalent is 

IxM/h(l) = [l - M) - "2(1 - M) 2 ] , 

2.2.2 LS2: least squares with 7(1) free 

Relaxing the constraint that 7a (1) = 7a. (1) gives laws that 
can be again be solved in a trivial least-squares exercise, 
with 7a (1) as an additional free parameter: 



7a(m) = 7\(1)[1-<1-|u)] 



(linear), 



7 A (At) = 7 A (1) [1 - ui(l - p) - u 2 (l - /i) 2 ] 

(quadratic; in practice, both sides may be divided by 7a (1) 
in these two equations). 

2.2.3 Flux-conserving fit: FC1 

The physical flux Fa is related to the specific intensity 
through 



Fa = 2tt f h(fi)fidfi 
Jo 

= 4tt77a 



(4) 



3 In practice, even 'two-parameter' fits still allow only one coef- 
ficient to be usefully constrained; see the discussion in Section]?] 



4 A number of other limb-darkening laws have been proposed; 
cf., e.g., Di'az-Cordoves & Gimenez ( 1992 i 
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Figure 1. Comparison of linear limb-darkening coefficients deter- 
mined from model atmospheres by different numerical techniques 
(cf. Section^, as a function of effective temperature. Open circles 
are LSI results. 



where H\ is the Eddington flux (the first-order moment of 
the radiation field). The integration of eqtn. fusing an an- 
alytical limb-darkening law to represent with coeffi- 
cients determined by least squares, will not normally recover 
the physical flux exactly. To address this, we can impose the 
condition that 



that is, 



Fx =2n I x (p)fJ.dfi- 



Fx = 



Fx = 



7ri(l) [1 - u/3] 

271-/(1) 



12 



6 - 2ui 



in the linear and quadratic cases, respectively. Requiring 
7a (/Li), evaluated from the limb-darkening law, to equal 
Ix(l-t-), evaluated from the model atmosphere, at some ar- 
bitrary fi = i, we obtain 



nlx{x) - Fx 



(tyIx(x)/3) + (x-1)F x 



for the linear law. |Wade fc Rucinski ( 1985 I chose x 
whence 

u = 3 [1 - F x / (tt/a(1))] 



(5) 
= 1, 



(noting the Wade & Rucinski s "angle-averaged" [astrophys- 
ical] flux is Fx/n in the nomenclature adopted here). In ef- 
fect, the choice of x fixes the intercept of the linear law, with 
the constraint of flux conservation then fixing the slope. 

The equivalent algebra for the quadratic law follows 
from selecting any two values /i = x\ , x^ at which Ix (/i) is 
equal to Ix(^), giving a pair of simultaneous equations that 
can readily be solved for Ui,U2- Wade & Rucinski (19851, 
and subsequent authors, used X\ = 1, X2 = 0.1 (values which 
are also adopted here), but again these are more or less ar- 
bitrary choices. 



2.2-4 Flux- conserving least squares: FC2 

The weakness of the standard flux-conserving approach is 
the lack of a compelling physical argument to select any 
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Figure 2. Limb darkening in the H and U bands for a 4 kK 
model. The thick 'lines' are individual model-atmosphere inten- 
sities, shown as points which merge together at this scale. The 
fitted 4-coefficient limb-darkening laws are shown drawn through 
the points. Straight lines show the linear limb-darkening laws for 
the H band, with coefficients determined by standard flux con- 
servation and by flux-conserving least squares (FC1, FC2, respec- 
tively; cf. Section|2.2|l. 



particular x values for the normalization (other than requir- 
ing the intensities to be everywhere positive; e.g., requiring 
^ u ^ 1 in the linear case). 

Rather than making an arbitrary choice of x, we can 
instead introduce the more objective requirement of min- 
imising the sum of the squares of the differences between 
model and fitted intensities while still requiring flux to be 
conserved. For a linear law it is convenient first to deter- 
mine u by by minimising (i (p) — I(p)j , using standard 
least-squares techniques, where 



Mm) = 



3Fx 



1 



- u(l - 
3 — it 



/') 



and to then evaluate 



3F X 



7r(3 — It) 

Corresponding results for the quadratic law are 



QF X 

7T 



• Ui(l — fl) — u 2 {l — fi) 2 



6 - 2ui 



■a 2 



h(l) = 



6Fx 

7r(6 — 2ui — U2) 

Not surprisingly, this newly introduced approach of flux- 
conserving least squares generally yields numerical coeffi- 
cients very close to those found using the LS2 method. 
Therefore, although it may be regarded as superior to LS2 
in principle, in practice it affords no great benefit (and turns 
out not to give results particularly close to photometrically 
inferred LDCs). 

2.3 Other numerical factors 

The foregoing numerical methods can (and do) yield sub- 
stantially different LDC values, even for the standard lin- 
ear and quadratic representations of a given, fixed, intensity 
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distribution, as is illustrated by Figs. [T] and [2] For a given 
intensity distribution, in the optical wavelength regime the 
FC1 u coefficient is usually the smallest numerically; LS2 
and FC2 u coefficients are very similar, and relatively large; 
and the LSI coefficient is intermediate. 

When characterizing model-atmosphere results, the 
density and distribution of angles at which intensities are 
calculated, and the weighting scheme, will also influence the 
numerical values of limb-darkening coefficients (e.g., |Dfaz- 



Cordoves fc Gimenez1|1992| |Claret|[2008| . In the present 
work, intensities were computed for ji = 0.001 to 1 at steps of 
0.00l(f] and equally weighted when fitting functional forms. 

Of course, the physics used in constructing the model 
atmosphere is also critical. The limb-darkening coefficients 
used throughout this paper were computed using the At- 
LAS9 line-blanketed LTE model-atmosp here code (|Kurucz 



1993), as ported to GNU-linux systems by Sbordone, Bonifa- 



cio & Castelli (20071, with the Opacity Distribution Func- 
tions described by |Howarth|p011[ ). Solar abundances, a mi- 
croturbulent velocity of Vt = 2 kms -1 , and mixing-length 
parameter 1/ H = 1.25 were adopted unless noted otherwise. 
These models use time- independent, plane-parallel struc- 
tures; while atmospheric extension is unlikely to be impor- 
tant in the parameter space discussed here, the neglect of 
time-dependent 3D effects may be significant when compar- 
ing with empirical results (e.g., Bigot et al.|2006 i. 



3 INFERENCES FROM EXOTRANSIT 
PHOTOMETRY 

The dispersion in coefficient values introduced simply by 
numerical techniques poses the question: which procedure is 
most appropriate for comparing model-atmosphere results 
with observational determinations of limb darkening? To an- 
swer this question it is necessary first to examine just what 
it is that is measured from transit observations. 

Photometric observations of exoplanetary transits 
record, essentially, the variation of 7a (r) / F\ along a chord. 
In practice, this variation is parameterized by an analytical 
limb-darkening law, whose coefficients are optimized as part 
of the global fitting procedure. This optimization process is 
quite different from fitting a limb-darkening law to model- 
atmosphere intensities, so it is immediately clear that there 
cannot be any simple one-to-one correspondence between 
photometric and model-atmopshere LDCs. 

Moreover, the extent to which a transit light-curve en- 
codes the global limb darkening must depend on the im- 
pact parameter]^] on how faithfully the chosen parametriza- 
tion of the limb darkening matches the intensities on those 
parts of the star that are occulted; and how well it extrap- 
olates to those parts that are not. Given that simple lin- 
ear and quadratic limb-darkening laws give only approx- 
imate representations of actual intensity distributions, it 
might be anticipated that LD coefficients determined from 
light-curves must, at some level, suffer systematic biasses 

5 The intensities were calculated in detail, not 'densified' from a 
sparser grid. 

6 The impact parameter is p = (a cos i)/R* for a circular orbit, 
where a is the orbital semi-major axis and i is the orbital incli- 
nation 



depending on impact parameter, reinforcing the point that 
these photometric coefficients must fail to correspond di- 
rectly to any equivalent, single- valued, characterization of 
model-atmosphere results. 

To demonstrate and quantify these effects, model tran- 
sit light-curves were first generated for a range of stellar tem- 
peratures and passbands. Limb darkening was represented 
in the light-curve calculations by eqtn. [3] with coefficients 
determined by least squares (i.e., for practical purposes, the 
model-atmosphere results were represented almost exactly). 
The geometry was set by choosing a ratio of planetary to 
stellar radii of 1:10 and a centres-of-mass separation of 10/?*, 
representative of 'hot Jupiter' systems commonly observed 
to transit (although the subsequent limb-darkening results 
are insensitive to the adopted values), with the impact pa- 
rameter varied over the range 0-0.95. 

The resulting light-curves were then solved for basic 
geometrical parameters, and for linear or quadratic LDCs. 
These calculations were performed using a modified ver- 

which is it- 



sion of jktebop (Southworth et al. 2004a 



self based on Etzel's ebop code ( Etzei|1981 



b) 



Popper & Etzel 



19811. Proximity effects (tidal distortion, 'reflection', etc.) 
were throughout assumed to be negligible. 

3.1 Linear law 

Figure [3] shows a selection of the results, and demonstrates 
that, in this parameter space, linear LDCs derived from pho- 
tometry systematically increase with increasing impact pa- 
rameter, by up to ~0.2 (~60%). This behaviour is straight- 
forward to understand: in the optical regime dl/dfi de- 
creases with increasing fj,, so that any transit that is not 
central (i.e., inclination i 7^ 90°) samples a relatively steep 
part of the limb-darkening law. A linear approximation to 
that law must therefore yield a linear limb-darkening coef- 
ficient that is, in general, larger than that derived from a 
central transit. 

Figure [2] illustrates this point by showing model- 
atmosphere results, and linear approximations, for a 4 kK 
model. The intensity in the U band is very nearly a linear 
function of /1, so any characterization will yield similar nu- 
merical values for the u coefficient. This is confirmed both 
in that least-squares and flux-conserving approaches yield 
similar results in this passband, and in that the u coefficient 
derived photometrically is insensitive to impact parameter 
(upper-left panel of Fig. |3j. This contrasts with //-band re- 
sults; the intensity there is a strongly non-linear function 
of angle, and any diagnostic that characterizes only small [i 
values must yield a larger u coefficient than that character- 
izing the entire centre-to-limb variation. Fig.[3]shows this to 
be the case. 

Fig. [3] also shows the linear limb-darkening coeffi- 
cients obtained by fitting the input intensity distributions 
directly, using flux-conserving (FC1) and flux-conserving 
least-squares (FC2) techniques, which bracket the range 
of numerical values derived directly from the model atmo- 
spheres. In this parameter space, these coefficients also al- 
most bracket the corresponding photometric LDCs which 
suggests a simple, if rough-and-ready, means of comparing 
observational and model-atmosphere parametrizations. Fur- 
thermore, if one had to choose a single, linear limb-darkening 
coefficient to compare with photometric results, then the 
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Figure 3. Photometrically determined linear limb-darkening coefficients (small dots), for fixed input limb darkening, in the Johnson- 
Cousins-Glass UVIHL passbands (cf. Section |3.1| . Larger symbols show corresponding linear limb-darkening coefficients determined 
directly from the same input model-atmosphere intensity distributions using flux-conserving and flux-conserving least-squares fitting 
(FC1, large squares, FC2, large diamonds; Section|2,2| These single-valued results are plotted at arbitrary impact parameters). 



FC1 value is probably the least poor option; for randomly 
inclined orbits, smaller impact parameters are more proba- 
ble than larger ones (and observational selection effects also 
favour higher orbital inclinations), and photometrically de- 
termined limb darkening coefficients are generally closest to 
the FC1 LDC in this caseQ 



onto new principal axes, 



11)1 = Ml COS (p — U2 sm ( 

w 2 — u 2 cos 4> + Ui sin < 



(6) 



3.2 Quadratic law 

The variations in linear limb-darkening coefficient are 
present a fortiori for the quadratic coefficients, although 
here the interpretation is less straightforward because of the 
well-known strong correlation between ui and u 2 , evident in 
the left-hand panels of Fig. |4| (see also Fig . 2 in Southworth 
20081). IPall (120081) and Kipping & Bakosl (|2011al) point out 



that this correlation is largely removed through a rotation 



7 This conclusion is supported by results from many more syn- 
thetic light-curves than are reported on here. 



with (f> ~ 40°. Results rotated to these co-ordinates are 
shown in the right-hand panel of Fig. [4j and confirm that, 
while the wi values continue to show a large variation with 
impact parameter, W2 is more nearly constant. 

The correspondence between the photometric and 
model-atmosphere results is less straightforward than with 
the linear law; the model-atmosphere representations of 
limb darkening show no simple relationship to the photo- 
metrically-determined equivalents (notwithstanding rough 
quantitative similarities). Nevertheless, the small dispersion 
found for the w 2 coefficient in both observational and model- 
atmosphere characterizations of limb darkening indicates 
that this should be the parameter of choice when making 
comparisons. 
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Figure 4. Photometrically determined quadratic limb-darkening coefficients for a 4 kK model in the Johnson-Cousins-Glass UVIHL 
passbands. Left-hand panels, Mi,M2 coefficients (eqtn.[2]|; right-hand panels, rotated w\,W2 coefficients (Section |3.2| . Larger dots show 
corresponding quadratic limb-darkening coefficients determined directly from the same input model-atmosphere intensity distributions 
using flux-conserving and flux-conserving least-squares fitting (FC1, large squares, FC2, large diamonds; Section |2.2| These single-valued 
results are plotted at arbitrary impact parameters). 



4 COMPARING MODEL-ATMOSPHERE AND 
PHOTOMETRIC LIMB-DARKENING 

The foregoing sections emphasize that different characteri- 
zations of model-atmosphere results can give quite different 
numerical results (e.g., Fig. [I]); and that light-curve analyses, 
using, of necessity, approximate limb-darkening 'laws', yield 
LDCs that vary with transit geometry (e.g., Fig. j3|). Fur- 
thermore, although a given analytical limb-darkening law is 
adopted in photometric studies, the determination of its co- 
efficients through light-curve modelling is, numerically, fun- 
damentally distinct from the techniques of fitting model- 
atmosphere intensity distributions discussed in Section [2] 
that is, comparing photometric and model-atmosphere re- 
sults is, to an extent, like comparing apples and oranges 
(but see Sandford 1995; Baronc 20001. In order to examine 
the relationship between empirical, photometric LDCs and 
theoretical model-atmosphere values, it is therefore neces- 
sary to devise a method ensuring a fair comparison. 

The most direct way to perform such a 'like for like' 
comparison is to adapt the methods used in Section [3j 
i.e., to generate model light-curves for well-studied systems, 
using as inputs the empirically determined geometric pa- 
rameters, coupled to model-atmosphere intensity distribu- 
tions (in practice, approximated by eqtn. |3| for the 'known' 



stellar parameters. This synthetic photometry can then be 
solved for the geometric parameters and LDCs, using the 
same simplified limb-darkening law adopted in the observa- 
tional photometric analysis. The resulting hybrid synthetic- 
photometry/atmosphere-model (SPAM) LDCs can reason- 
ably be compared directly with empirical values]^] 

This approach has been used to investigate two illustra- 
tive datasets: the eight stars with Kepler data analysed by 
(|2011a|b|), and the multiwavelength study 
Knutson et al. (2007 see also |Southworth| 



Kipping & Bakos 



of HD 209458 by 



2008 Claret 2009). Synthetic light-curves were generated 



with 1000 data points through transit (phases ±0.05). Sta- 
tistical errors are not quoted on any results because the anal- 
ysis is essentially deterministic. 



4.1 Kepler targets 



|Kipping fc Bakos| ( 2011a|b l derived quadratic limb- 
darkening coefficients for the eight Kepler targets they stud- 



8 Using a simplified description of limb darkening in the fitting 
step drives the inferred geometric parameters away from their 
input values, but by usually unimportant amounts (AppendixlBII. 
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frames because this figure is a projection of m ultiparameter modelling ont o a specific 2D plane). Green squares show the median values 
from MCMC runs (the solutions adopted by |Kipping fc Bakos|[2011a b|, and green diamonds the minimum-x 2 MCMC results. Red 



dots show fits to model-atmosphere intensities (left to right: FC2/LS2 [indistinguishable at this scale], LSI, FC1), while horizontal and 
vertical lines indicate the SPAM solutions. The small rectangle shown in each panel (perhaps most easily seen by zooming in on the 
on-line version) encompasses the SPAM solutions for all six targets, and is included to provide an qualitative indication of the rather 
small scale of uncertainties likely to result from any plausible errors in input stellar parameters. The rotated w\,W2 axes (eqtn. [6| are 
shown in the TrES-2 panel, for reference; by design, most of the variance in the MCMC results is in w±. 



ied. Their Monte-Carlo Markov-Chain results are repro- 
duced here in Fig. [5] 



Custom model atmospheres were computed for each 
system as described in Sectio n |2.3| adopted stellar parame- 
ters are summarized in Table [Al| (Appendix |A|). This group 
of stars samples a fairly small range in atmospheric prop- 
erties (T off = 5647:6297 K, logg = 3.96:4.59, [M/H] = 
— 0.55:+0.33), which is reflected in a rather small range in 
model-atmosphere and SPAM LDCs. It's therefore some- 
what surprising that agreement between empirical LDCs 
and those from the SPAM approach (or from direct fitting 



to model atmosphere^ varies from excellent (Kepler-6) to 
statistically unacceptable (e.g., Kepler-5). 

There is a suggestion that the extent of agreement cor- 
relates with temperature; the SPAM LDCs for the three 
coolest stars fall within the cloud of the best-fitting 90% 
of solutions, while those for the three hottest lie (just) out- 
side [^] The trend is for the cloud of empirical values to move 
towards smaller (ui, U2) values with increasing temperature, 



9 In general, the SPAM results are closest to the FC1 direct char- 
acterization of intensities. 

10 The cooler stars in this sample are also those with higher grav- 
ities and metallicities, so temperature is not necessarily the key 
parameter. 
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Figure 6. Comparison of linear limb-darkening coefficients deter- 
mined observationally for HD 209458 (lower set of black points, 
from |Southw orth 2008) and SPAM calculations (upper set of red 
points). The 'error bars' on the model-atmosphere results are the 
result of varying input stellar parameters (see Section |4.2| for de- 
tails); for both SPAM and empirical results, the error bars are 
smaller than the points at most wavelengths. Continuous and 
dotted lines connect LS2 and FC1 model-atmosphere results, re- 
spectively. 
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Figure 7. Comparison of quadratic limb-darkening coefficients 
for HD 209 158. Black points with error bars, empirical values from 
|Southworth| l |2008| ); red points, SPAM results (Section [P). Con- 
tinuous and dotted lines connect LS2 and FC1 model-atmosphere 
results, respectively. 



compared to the model-atmosphere results. It's unclear why 
the empirical results should show so much greater variation 
than the models, suggesting that this apparent trend may 
simply be an artefact of the small sample, or that some ad- 
ditional factor plays an unexpectedly important role. 

4.2 HD 209458 

Baseline parameters of T e g = 6113 K ([Casagrande et al.| 



M/H] = +0.03 (Sousa et al. 20081, 



2010], log g = 4.50, 
Vt = 2 km s _ , i/H = 1.25 were adopted to construct the 
reference model atmosphere and intensities for HD 209458. 
Broad-band limb-darkening was calculated by assuming 'top 



hat' response functions for the photometric passbands of the 
|Knutson et al.| HST observations. The principal results are 
summarized in Table lA2l 

Additional models were run for T e g = 5913, 6313; 
log 5 = 4.2,4.8; t/H = 0.5; v t = 0,4 kms" 1 ; and [M/H] = 
—0.4, +0.4. These ranges allow for quite generous uncertain- 
ties in parameters for this well-studied system. The extremes 
in linear LDCs from the models are for the low-Teff and high- 
gravity models (numerically largest and smallest coefficients, 
respectively), and these models are used to illustrate plau- 
sible 'error bars' on the SPAM coefficients in Figs. [6] and [7] 

Fig. [6] shows results for linear coefficients. The discrep- 
ancies between model-atmosphere and photometric results 



already noted by Claret (20091 see also Southworth 20081, 



on the basis of older models persist in the new analysis. 

The comparison for quadratic coefficients is shown in 
Fig. [7] The variation with wavelength of both U\ and ui 
coefficients is much less for the SPAM coefficients than is 
found empirically. However, both sequences run almost par- 
allel to the rotated w\ axis, and agreement in the better- 
determined W2 parameter is tolerable at all wavelengths. [^] 
In particular, for the ~678 nm passband, which is close to 
the effective wavelength of the Kepler results, the agreement 
is reasonably good, [(w 1 ,w 2 ) = (0.234, 0.385), (0.099, 0.363) 
for SPAM and light-curve coefficients, respectively]. This is 
in contrast to the Kepler results for stars at similar effective 
temperatures (but is consistent with the result that it is the 
higher-gravity stars that show the best agreement between 
models and observations). 



5 SUMMARY AND CONCLUSIONS 

Different methods of fitting a given limb-darkening law 
to a given model-atmosphere intensity distribution lead to 
quite different numerical coefficients. Furthermore, the limb- 
darkening coefficients determined from photometry of exo- 
planctary transits are functions of impact parameter, and 
can't reliably be compared directly to any of the standard 
model-atmosphere characterizations. 

A more direct comparison can be made if the model 
intensities are translated into observer space, through the 
medium of synthetic light-curves. The resulting synthetic- 
photometry/atmosphere-model (SPAM) limb-darkening co- 
efficients are not single- valued, but can be compared directly 
with empirical results. 

[If one had to choose a traditional single- valued repre- 
sentation of model-atmosphere results, then at optical wave- 
lengths closest agreement with the SPAM results is gen- 
erally obtained with the standard (FC1) flux-conserving 

11 The principal cause of the minor quantitative differences be- 
tween the results shown in Fig. [6] appears to be the different 
treatments of convection adopted here and by |Kurucz| ( |1993| the 
source of Claret's models). 

12 The observational uncertainties can't be straightforwardly 
propagated from a{u\^) to o(w\ } 2) because of the strong corre- 
lation between u\ and 112 ■ It should also be noted that the error 
bars in Fig. [7] represent the 68% dispersion in parameter values 
obtained from MC replications (Southworth, personal communi- 
cation); they therefore can't be compared directly to the disper- 
sion in results shown in Fig. [5] which show parameters from the 
90/95% of solutions with the smallest overall x 2 values. 



Limb darkening and exotransits 9 



method, which also yields the smallest value of the linear 
limb-darkening coefficient.] 

For the commonly used quadratic limb-darkening law, 
most of the variation in different fits to model-atmosphere 
intensities is in the Wi parameter, with much smaller disper- 
sions in the W2 coefficient. Since the wi axis is also defined 
as that which maximizes dispersion in observational (Monte- 
Carlo) results, the most sensitive comparison between mod- 
els and observations is in W2- 

New model-atmosphere calculations, analysed with the 
SPAM approach, show mixed results. Agreement with em- 
pirical Kepler LDCs is good in some cases (differences in W2 
less than 0.06 in four out of six systems), but not in oth- 
ers. There is a hint of a possible temperature dependence 
in the extent of disagreement for these targets, with cooler 
stars showing better agreement. However, at similar effective 
wavelengths HST results for HD 209458 (which is at the hot- 
ter end of the range of Kepler targets) agree well with mod- 
els; there are discrepancies at longer and short wavelengths, 
though again with fair agreement in W2- Since gravity (and 
metallicity) correlate with temperature for the Kepler sam- 
ple, and since HD 209458 is both high-temperature and high- 
gravity in the context of that sample, this might be taken 
as an indication that agreement is better at higher gravi- 
ties (with temperature as a secondary factor) . However, the 
ranges in all quantitities characterizing the atmospheres are 
so small as to render such conclusions speculative at this 
stage. Forthcoming results from Kepler, and other missions, 
should enlarge the parameter space, and permit better dis- 
crimination of where models and observations do and do not 
agree. 
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APPENDIX A: FIT RESULTS 

APPENDIX B: SYSTEM ATICS OF 
GEOMETRIC PARAMETERS 

Because light-curves have only a rather weak dependence on 
limb darkening, we might expect that the use of simplified 
limb-darkening laws, or even moderately inaccurate LDCs, 
should have only a very modest effect on the determination 
of basic geometric parameters when modelling photometry. 
To demonstrate this (at the referee's suggestion), additional 
results from the model grids described in Section [3] are pre- 
sented in Fig |B1| 

[To remind the reader, model light-curves were gener- 
ated using a ratio of planetary to stellar radii of 1:10 and 
a centres-of-mass separation of 10/?*, over a range of im- 
pact parameters, and an effectively exact representation of 
limb darkening. These light-curves were then solved, adopt- 
ing linear or quadratic limb darkening (with the LDCs as 
free parameters). It is the results of these light-curve solu- 



tions that are summarized in Fig Bl ] 

The systematic errors in fitted geometric parameters 
(up to ~4% in r/R t and (r + R r )/a for the models presented 
here) might be significant for the best-quality photometry if 
a linear limb-darkening law is assumed, but are negligible if 
a quadratic law is used. 

This paper has been typeset from a T^X/ DTjtX file prepared 
by the author. 
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Table Al. Limb-darkening results: Kepler photometry. For each star, the physical parameters used in the model-atmosphere calculations 
are first listed, followed by the resulting 4-parameter Kepler-band limb-darkening coefficients (eqtn. [3J. Subsequent columns list the 
SPAM coefficients, determined by fitting synthetic light-curves generated from the 'known' system parameters; the empirical photo- 
metric coefficients; and, for reference, coefficients determined directly from the model-atmosphere intensity distributions. All necessary 
stellar & system parameters are adopted from |Kipping fc Bakos] l|2011a|b|| . 
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Table A2. Limb-darkening results: HD 209458. System parameters and observed LDCs are from |Southw orth (2008}. Stellar parameters 
T eS = 6113 K jCasagrande et al.|20l6"} , logg = 4.50, [M/H] = +0.03 jSousa et aL|2008) , v t = 2 krns" 1 , l/H = 1.25. Columns follow 
the model of Table |A1[ with rows grouped by HST wavelength in nm. 
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Figure Bl. Geometric parameters determined from modelling synthetic light-curves, as functions of input impact parameter. The 
light-curves were generated using essentially exact representations of limb darkening; ratio of planetary to stellar radii r/R- t = 0.1; the 
sum of the radii, in units of the semi-major axis r + R„)/a = 0.11. Light-curve solutions assume linear (left-hand panels) or quadratic 
(right-hand panels) limb darkening, with coefficients optimised as part of the fitting process. The (generally small) departures from input 
parameters are solely a result of using these approximate representations of limb darkening. Results are shown for stellar models at 4, 
6, 8, and 12 kK, and UBVRIJHKL passbands; 4-kK C/-band and 8-kK B-band results bracket most of the data when a linear law is 
adopted, and are highlighted with larger symbols. Results for the R band, which are representative of typical unfiltered CCD systems, 
are shown in red. 



